The Energetics of Li Off-Centering in Ki-^^Li^^TaOa: First Principles Calculations 
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Ki-iLiajTaOs (KLT) solid solutions exhibit a variety of interesting physical phenomena related to 
large displacements of Li-ions from ideal perovskite A-site positions. First-principles calculations for 
KLT supercells were used to investigate these phenomena. Lattice dynamics calculations for KLT 
exhibit a Li off-centering instability. The energetics of Li-displacements for isolated Li-ions and 
for Li-Li pairs up to 4th neighbors were calculated. Interactions between nearest neighbor Li-ions, 
in a Li-Li pair, strongly favor ferroelectric alignment along the pair axis. Such Li-Li pairs can be 
considered "seeds" for polar nanoclusters in KLT. Electrostriction, local oxygen relaxation, coupling 
to the KT soft-mode, and interactions with neighboring Li-ions all enhance the polarization from Li 
off-centering. Calculated hopping barriers for isolated Li-ions and for nearest neighbor Li-Li pairs 
are in good agreement with Arrhenius fits to experimental dielectric data. 



I. INTRODUCTION 

Potassium tantalate, KTaOs (KT), is a cubic per- 
ovskite with a quantum paraelectric ground state. ^ In 
Ki-ajLia^TaOs (KLT) solid solutions, Li-ions substitute 
for K on perovskite A-sites. Li-ions are smaller than 
K, and displace from ideal, centrosymmetric, A-sites by 
about 1.26 A along cubic [0,0,l]-type vectors ([0, 0, l]c)^ 
Displacements of Li-ions generate strong local dipole mo- 
ments [Li-dipole(s)] which couple electrostatically to the 
KT polar soft-mode. Near neighbor (nn) Li-dipoles can 
also interact to form polar nanoclusters (PNC); i.e. mul- 
ticell regions in which the local polarization direction is 
strongly correlated. The complex nature of ferroelectric 
(FE) ordering in KLT derives from a combination of in- 
teractions between Li-dipoles, the soft-mode, and PNC. 

In the composition range < x < 0.02, KLT is a quan- 
tum paraelectric, but for 0.02 < x < 0.06, both FE and 
dipole-glass characteristics are observedi^i^ In the tem- 
perature (T) interval < T < 100 K, the dielectric per- 
mittivity, e{T), has a peak, and second harmonic gener- 
ation exhibits an abrupt increase, with hysteresisAi The 
existence of randomly distributed and randomly poled 
PNC was established from dielectric measurements, bire- 
fringence and nuclear magnetic resonance (NMR)i^i^*24S 
Neighboring PNC may touch, but PNC-coarsening is 
apparently prevented by random long-range PNC-PNC 
interactionsii In zero applied field, a glassy phase is 
observedf^ however, a sufficient applied field will induce 
FE long-range order (LRO).^ For 0.06 < a; < 0.15, 
KLT undergoes a FE phase transition at T > lOOK, even 
in the absence of an applied field. At x > 0.15, the per- 
ovskite solid solution is no longer stable. 

Dielectric measurements on KLT with < a; < 
0.02 exhibit only one relaxational peak in e{T), while at 
X « 0.04, two peaks are observed, and they are associated 
with two relaxational processes (fast and slow)i^ii24iiiiS 
Temperature variations of both relaxation processes fol- 
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1000 K (86 meV), and the slow pro- 



cess, at 0.02 < X < 0.06, has been reported as 2100 K^ 
2400 Kfii and 2800 K.^^ 

The fast process 1000 K activation energy is believed 
to be the barrier for an isolated Li-ion to hop from one 
[0,0, l]c position to another separated by a 90° angle. 
NMR data support this hypothesis^ The slow process 
with an activation energy of 2100 <U< 2800 KiitiSJ^ 
is less well understood. On the basis of acoustic mea- 
surements it was hypothesized®'^" that the slow process 
might be 180° correlated Li-Li nn pair reorientations (e.g. 
z — z ^ z — 'z relaxation for nn separated by aoz where 
ao is the lattice parameter of the primitive cubic unit 
cell, and £ is a unit vector in the [0, 0, z] direction). 

The energy spectrum and kinetics of PNC reorienta- 
tion is a general problem for FE systems that exhibit 
both relaxor and soft-mode behavior (cfii^). KLT is 
a relatively simple model system for such phenomena 
because isolated Li-ions in KLT solutions can be re- 
garded as randomly distributed sources of 6-state Potts- 
like dipole fields piStiSiii Previous modeling of KLT in- 
cludes: Semi-empirical shell model calculationaiS; In- 
termediate neglect of the differential overlap (INDO) 
calculations^^'^*'; Full potential linear muffin-tin orbital 
(FPLMTO) calculationE^°-^^ (without ionic relaxations). 
All these studies yielded potential barriers for 90° Li- 
hopping that were significantly lower than those deduced 
from experimental fits. 

KLT solutions with 0.016 < a: < 0.05 exhibit a 
large photocurrent below 80Ki ^^i^'^i^^ No photocurrent 
has been observed in pure KT, nor in solid solutions of 
KT with KNbOs. Photocurrent is usually thought to 
be associated with 0-hole centers near Li-occupied A- 
sites^^ . Recent INDO calculationsSS confirmed that dis- 
placements of 0^^-ions which are nn of Li create shallow 
states in the forbidden gap. 

First-principles (FP) computations allow determina- 
tions of the electronic structures and energetics as func- 
tions of ionic coordinates, and thus provide a sound basis 
for understanding and quantitatively modeling KLT. In 
this work, FP density functional theory (DFT) methods 
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were used to calculate the energetics of Li off-centering 
for isolated Li-ions, and for Li-Li pairs up to fourth neigh- 
bors. The energetics of 90° isolated Li-ion hopping and 
180° Li-Li nn-pair hopping were calculated. Lattice dy- 
namics studies of KLT demonstrate correlated motions of 
Li and the surrounding ions (particularly nn O^^). Total 
polarization from Li off-centering, and individual contri- 
butions from different ions were quantified. Electronic 
band structures were calculated to investigate shallow 
electronic states that may promote photocurrent. 



II. METHODS 

All DFT calculations were done with the Vienna ab 
initio simulation package (VASP)2S*2&. A plane wave ba- 
sis set for electronic wavefunctions and ultrasoft pseu- 
dopotentials were used^^ in the local density approxima- 
tion (LDA) for exchange and correlation energies. VASP 
computes interatomic forces and total energies for crys- 
tals, and allows either global or constrained relaxations of 
internal coordinates and/or lattice parameters. Frozen- 
phonon methods were used to obtain force constants for 
computing lattice dynamics. Berry phase analyses, as 
implemented in VASP by M. Marsman, were used to cal- 
culate dynamical charges. 

Investigating the effects of Li-ions in KLT with low Li 
concentration requires large supercells, e.g. (Fig. 1): (a) 
KT40, a 40 atom supercell of pure KTaOs; (b) KLT40, 
a 40 atom supercell with basis vectors [0,0, 2]c, [0, 2,0]c 
and [2,0,0]c (units of uq, x = 0.125); (c) KLTSOa an 
80 atom supercell with basis vectors [2,2,0]c, [2,0, 2]c, 
and [0,2,2]c and x = 0.0625; (d) KLT806, same as 
KLT80a but with two Li per supercell, x = 0.125, sep- 
arated by aoz; (e, f, g) KLT80c, KLT80fi, KLTSOe, same 
as KLT80b, but with two Li separated by aQ{x + y), 
ao{x + y + z) and 2aoz, respectively. For KT40 and 
KLT40, a 4x4x4 Monkhorst-Pack fc-point grid was used 
and for KLTSOa, etc., a 2x2x2 grid. Results for KLT80a 
are negligibly different with a 4x4x4 Monkhorst-Pack 
grid. 

This computational scheme was tested on pure KT 
with 5 atoms per cell. The VASP lattice constant a 
= 3.96 A agrees with previous LDA computations^*, 
and is smaller than the experimental values (between 
3.987 A and 3.989 A at room temperaturo^i^?!??!?^). 
The calculated bulk modulus for a KT40 supercell with 
4x4x4 fc-point grid is 2.0 Mbar, which is lower than the 
experimental^ value 2.3 Mbar. A first principles LDA 
result for a 5-ion primitive cell with an 8x8x8 fc-point 
grid is 2.1 Mbar—, but a full-potential linear muffin-tin 
orbital calculation, primitive cell and 8x8x8 fc-point grid 
yieldedji^ 2.3 Mbar. Hence, ultrasoft pseudopentials pro- 
vide reasonable accuracy compared with other ab initio 
calculations. 



III. RESULTS 
A. Reference Structures 

To study Li off-centering, one needs a reference struc- 
ture in which Li occupies an ideal centrosymmetric A- 
site. A 40-atom reference structure, KLT40Rof, is ob- 
tained by placing all ions of KLT40 on ideal perovskite 
positions and optimizing lattice parameters and inter- 
nal coordinates (except Li), subject to a cubic symmetry 
constraint. The resulting cell has a = 2 x 3.956 A. Crys- 
tallographic data for KLT40Rcf , are tabulated in TableQ] 
Because Li is smaller than K, there is a volume contrac- 
tion around Li such that nn Ta ions displace by 0.010 
A towards Li, and nn O ions displace by 0.024 A to- 
wards Li. 

TABLE L Calculated structure for KLT40fle/:KLT40 with Li 
occupying ideal A-site positions, and the structure relaxed, 
subject to a cubic symmetry constraint. Space group PmSm; 
a = 2 X 3.9559 A. 



ion 


WyckofT position 


X 


y 


z 


Ki 


la 


0.0000 


0.0000 


0.0000 


K2 


3d 


0.0000 


0.0000 


0.5000 


Kg 


3c 


0.5000 


0.5000 


0.0000 


Li 


lb 


0.5000 


0.5000 


0.5000 


Ta 


8g 


0.2508 


0.2508 


0.2508 


Oi 


12i 


0.2509 


0.2509 


0.0000 


O2 


12j 


0.2522 


0.2522 


0.5000 



A series of KLT40R,ef frozen phonon calculations, with 
each symmetry-independent ion displaced in turn, was 
used to compute force constants, construct the dynamical 
matrix, and compute zone-center normal mode frequen- 
cies and eigenvectors. Symmetry analysis of KLT40Rof 
zone-center phonons yields 3Aig + 2A2g + 3A2„ + bEg + 
3Eu + 5Fig + UFiu + + 7F2u- Normal mode fre- 
quencies are listed in Table^J One zone-center instabil- 
ity is found for KLT40R,of, with symmetry Fi^. Thus, as 
expected, the high-symmetry reference structure is dy- 
namically unstable. Table IIIII gives the Fiu instabil- 
ity eigenvector, which is dominated by a Li off-centering 
[0 , , d] c displacement . 

For comparison, the force constants of a 40-atom cell of 
pure KTaOs (KT40) were calculated at the same lattice 
parameter. The TO normal mode results are listed in 
Table |TV| Force constants for KT40 and KLT40Ref are 
nearly identical. The only interatomic force constants 
that change by more than 0.41 eV/A^ are those involv- 
ing the 12 0-ions that are nn of Li. For example, the 
radial interionic force constant between an A-ion and a 
nn-0 changes sign from -0.21 eV/A^ in KT to -1-0.74 
eV/A^ in KLT40Ref when K is replaced by Li. The posi- 
tive sign for this term in KLT40Rcf implies that Li-motion 
in opposition to its nn 0-ions is energetically favorable. 
Thus, the single lattice instability in KLT40Rcf is domi- 
nated by Li-motion (85% of the dynamical matrix eigen- 



3 



TABLE II: Computed zone-center normal mode symmetries 
and frequencies {f in cm~^) for KLT40Rot (TableHJ. Fiu fre- 
quencies are for transverse optical modes, except for the zero- 
frequency acoustic triplet. Also listed are relative amplitudes 
arei for each mode in the fully relaxed tetragonal KLT40 
structure. The relaxation is dominated by the eigenvector 
of the 172 i cm~^ instability. In addition to the ferroelectric 
distortion itself which has Fiu symmetry, there are secondary 
distortions of Aig- and _Eg-type. 



symm. 


V 




symm. 


V 




symm. 


V 




Alg 


211 


0.0906 


Fig 


70 


0.0000 


F2g 


63 


0.0000 


Alcj 


418 


0.0877 


Fig 


205 


0.0000 


F2g 


178 


0.0000 


Alg 


445 


0.0504 


Fig 


338 


0.0000 


F2g 


262 


0.0000 








Fig 


527 


0.0000 


F2g 


327 0.0000 


Mg 


167 


0.0000 


Fig 


611 


0.0000 


F2g 


534 


0.0000 


A2g 


285 


0.0000 








F2g 


910 


0.0000 








Flu 


172 i 


2.9599 








A2u 


58 


0.0000 


Flu 





0.0000 


F2u 


129 


0.0000 


A2u 


507 


0.0000 


Flu 


120 


0.3047 


F2u 


159 


0.0000 


A2u 


945 


0.0000 


Flu 


164 


0.0596 


F2u 


178 


0.0000 








Flu 


177 


0.0080 


F2u 


243 


0.0000 


E, 


125 


0.6230 


Flu 


187 


0.2049 


F2u 


256 


0.0000 


E, 


213 


0.1954 


Flu 


200 


0.0082 


F2u 


331 


0.0000 


Eo 


279 


0.2125 


Flu 


205 


0.2038 


F2u 


349 


0.0000 


E, 


422 


0.1413 


Flu 


244 


0.0961 








E, 


448 


0.1241 


Flu 


336 


0.0901 














Flu 


364 


0.0041 








Eu 


69 


0.0000 


Flu 


432 


0.0373 








Eu 


495 


0.0000 


Flu 


555 


0.0378 








Eu 


577 


0.0000 


Flu 


862 


0.0075 









TABLE III: Dynamical matrix eigenvector for the z-polarized 
172 i cm~^ mode. Ionic labels correspond to those in tabled] 



ion 


x 


y 


z 




ey 




Ki 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0206 


K2 


0.0000 


0.0000 


5.0000 


0.0000 


0.0000 


-0.0374 


K2 


0.5000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0036 


Kg 


0.5000 


0.5000 


0.0000 


0.0000 


0.0000 


0.0507 


Ks 


0.0000 


0.5000 


0.5000 


0.0000 


0.0000 


-0.0350 


Li 


0.5000 


0.5000 


0.5000 


0.0000 


0.0000 


0.9252 


Ta 


0.2508 


0.2508 


0.2508 


0.0300 


0.0300 


0.0252 


Ol 


0.2509 


0.2509 


0.0000 


0.0000 


0.0000 


-0.0509 


Ol 


0.0000 


0.2509 


0.2509 


-0.0187 


0.0000 


-0.0301 


O2 


0.2522 


0.2522 


0.5000 


0.0000 


0.0000 


-0.0242 


O2 


0.5000 


0.2522 


0.2522 


0.0000 


-0.0649 


-0.0875 



vector) opposite its nn O (12% of the dynamical matrix 
eigenvector). Because KLT40Rof and KT40 have such 
similar force constant matrices, their normal mode spec- 
tra are very similar. Note that each TO mode in KT40 
has a corresponding TO mode of very similar frequency 
in KLT40Ref.^ 



TABLE IV: Computed TO normal mode frequencies {u in 
cm~^) and dynamical matrix eigenvectors for KTaOa at a = 
3.956 A. 



f 


UK 


UTa 


uo\\ 


uo± 


115 


-0.2989 


0.5419 


-0.3937 


-0.4807 


205 


0.8745 


-0.1737 


-0.2401 


-0.2715 


555 


0.0012 


0.0341 


-0.8531 


0.3682 



B. Fully Relaxed Structures 

Full ionic relaxation in KLT40 indicates that (consis- 
tent with previous unrelaxed computations^^^ii) maxi- 
mum energy reduction occurs when Li-ions are displaced 
to [0,0, (i]c (Figs. 2-3). The calculated equilibrium Li- 
displacemenfe^ is ci=1.009 A (Table EJ, whereas fits to 
experimental data yield^ d = 1.26 A. For comparison, 
the curves obtained when all ions, except for Li, are fixed 
at ideal perovskite positions, are also shown in Fig. 3. 
Clearly, the depth of the well and the magnitude of Li- 
displacement are sensitive to relaxations of the surround- 
ing ions, consistent with earlier results^ 

The largest sympathetic distortion associated with 
Li-displacements involve four nn 0-ions (see Table Ivl 
and Fig. 2), which are displaced towards Li by 0.109 
A. Similar correlated Li- and 0-displacements occur in 
LiNbOs and LiTaOs^^^. In each Li displaces, 

the oxygen atoms that Li approaches move toward the 
Li, and those left behind back away, as might be ex- 
pected from electrostatics. The magnitude of Li dis- 
placement is about 0.3 A in LiTaOs — or 0.4 A with 
full geometry optimization,-3& i.e., much smaller than Li 
in KLT. Note, however, that LiTaOa and LiNbOs have 
the LiNbOs structure type, a derivative of the corundum 
structure, not the perovskite type, and that the displace- 
ments patterns are therefore not directly comparable. 

The results listed in Table IVI indicate that Ta-ions are 
generally displaced in the same direction as Li, while O- 
ions shift in the opposite direction. This is consistent 
with the view that Li-displacement leads to freezing of a 
soft-mode polarization fluctuation that has a large corre- 
lation radius. 

A breakdown of the displacement pattern into normal 
mode coordinates for KLT40R,of is shown in Table [HJ 
Anharmonic coupling leads to some contributions from 
higher frequency modes, but otherwise, the displacement 
pattern is dominated by the eigenvector of the struc- 
tural instability. Vibrational frequencies for fully relaxed 
KLT40 were also calculated, and symmetry analysis of 
the fully relaxed tetragonal structure indicates that the 
phonon spectrum is 22 Ai + 8A2 + 14_Bi + 12^2 -I- 32^;. 
Modes belonging to irreducible representations Ai and 
E are infrared active (except for the Ai + E acoustic 
set). The frequencies of infrared active modes are listed 
in Tables IVll Each column represents the (summed) pro- 
jection (in percent) of the dynamical matrix eigenvec- 
tor of relaxed KLT40 onto one (or more) eigenvectors of 
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KLT40Rcf • If two eigenvectors were identical, the entry 
would be 100. 

The KT soft-mode frequency (at equilibrium lattice 
constant), 120 cm^^, splits into two modes in fully re- 
laxed KLT40: E (124 cm'^), and Ax (163 cm"!). 
Hence, the averaged frequency for these modes is in- 
creased by Li-displacement as indicated by experimental 
data for hardening of the soft-mode in KLT;*^-. However, 
calculated hardening of the soft mode is only significant 
for the component of the KT soft-mode triplet that is 
polarized in the Li-displacement direction. 

To study x-dependence of the soft mode splitting, lat- 
tice dynamics were calculated for KT40 (x = 0), KLTSOa 
{x = 0.0625), and KLT40 {x = 0.125), aU with fixed 
a — 3.9879A (the experimental value'^^ for x = 0.05). 
Results for the soft mode splitting are shown in Ta- 
ble The calculated splitting increases with x. TO 
splitting in the ferroelectric phase of KLT has been ob- 
served experimentally by Raman spectroscopy'^® and by 
neutron scattering^. The measured splitting'^^ in tetrag- 
onal KLT with x = 0.05 is given in Table EIIl The ex- 
perimental splitting is less than that obtained from FP 
results, owing to the use of periodic boundary conditions 
in FP calculations, which fixes the arrangement of the Li- 
sites and artificially aligns Li-dipoles from one supercell 
to the next (cf. Ref. f^]). 

The Li-dominated instability of KLT40R,of mixes with 
several other modes; its weighted average frequency is 
about 260 cm~^ in the polar direction and 330 cm^^ in 
transverse directions. 

To investigate the Li-isotope effect in KLT, phonon 
frequencies in Table IVll were recalculated using ^Li (to = 
6.015 amu) instead of natural Li (to — 6.941 amu). 
Individual mode frequencies increase when ^Li is used. 



TABLE V: Computed fully relaxed structure of KLT40. 
Space group P4mm; a = 2 x 3.9565 A; c = 2 x 3.9638 A; 
c/a ~ 1.002. Deviations are relative to positions in KLT40Rcf . 
The origin is chosen so that both structures have the same 
center of mass. WyckofI positions (pos.) are used. 



ion 


pos. 




X 




y 




z 


Sx, A 5y, A 


5z, A 


Ki 


la 





0000 





0000 





0014 


0.000 


0.000 


0.011 


Ka 


2c 





5000 





0000 





9996 


0.000 


0.000 


-0.003 


K3 


2c 





5000 





0000 





4963 


0.000 


0.000 


-0.029 


K4 


la 





0000 





0000 





4952 


0.000 


0.000 


-0.038 


Ks 


lb 





5000 





5000 





0074 


0.000 


0.000 


0.058 


Li 


lb 





5000 





5000 





6273 


0.000 


0.000 


1.009 


Tai 


4d 





2513 





2513 





2526 


0.004 


0.004 


0.015 


Taa 


4d 





2494 





2494 





7499 


-0.011 


-0.011 


0.005 


Oi 


4d 





2515 





2515 





9935 


0.005 


0.005 


-0.051 


O2 


4d 





2474 





2474 





4957 


-0.038 


-0.038 


-0.034 


O3 


4e 





0000 





2500 





2436 


0.000 


-0.007 


-0.058 


O4 


4e 





0000 





2533 





7457 


0.000 


0.019 


-0.027 


O5 


4f 





5000 





2482 





2442 


0.000 


-0.031 


-0.063 


Oe 


4f 





5000 





2589 





7379 


0.000 


0.053 


-0.079 



TABLE VL Computed TO normal mode frequencies (y in 
cm~^) for modes of symmetry E (left) and Ai (right) in 
fully relaxed KLT40. Also shown is the squared projection (in 
percent) of the eigenvector of each mode onto an eigenvector, 
or set of eigenvectors, of the dynamical matrix for KLT40Rcf. 
(1) The 172i cm"^ instability, (2) the 120 cm mode similar 
to the KT40 soft-mode. (3) All other polar eigenvectors and 
(4) All nonpolar eigenvectors. 



V 


(1) 


(2) 


(3) 


(4) 




(1) 


(2) 


(3) 


(4) 


66 


0.3 


0.6 


0.1 


99.0 


147 


0.1 


0.4 


17.2 


82.4 


79 


0.1 


0.1 


0.1 


99.7 


163 


0.1 


87.9 


8.4 


3.6 


124 


7.2 


89.0 


1.2 


2.7 


169 


0.8 


3 


1 


84.1 


12.0 


147 


0.0 


0.9 


0.8 


98.2 


176 


9.9 





1 

i 


85.6 


4.3 


159 


0.0 


0.2 


21.7 


78.1 


179 


9.4 





1 

i 


84.2 


6.3 


160 


0.6 


0.6 


77.6 


21.3 


199 


4.4 





1 


94.2 


1.2 


176 


0.1 


0.4 


23.2 


76.3 


205 


7.0 






1 

i 


59.9 


30.0 


178 


0.0 


0.1 


80.5 


19.4 


210 


2.8 


1 
i 


r 



14.8 


80.8 


187 


1.3 


0.9 


91.4 


6.4 


215 


0.2 


1 
i 


9 


35.5 


62.5 


199 


0.5 


0.6 


91.6 


7.3 


235 


22.4 








60.7 


16.6 


199 


0.5 


0.3 


51.7 


47.5 


262 


30.8 


1 
i 





57.5 


10.6 


206 


3.3 


1.1 


46.6 


49.1 


290 


5.0 





2 


2.2 


92.7 


224 


0.2 


0.1 


10.2 


89.5 


337 


2.2 





1 


96.7 


1.1 


239 


7.3 


0.6 


70.7 


21.4 


362 


0.0 








96.7 


3.3 


247 


1.6 


0.1 


6.2 


92.2 


417 


0.0 








24.4 


75.6 


252 


6.1 


0.3 


20.8 


72.8 


428 


2.8 








0.9 


96.3 


262 


0.0 


0.0 


0.7 


99.3 


444 


0.4 








63.9 


35.7 


307 


55.1 


3.4 


17.7 


23.9 


449 


1.4 








0.8 


97.7 


332 


0.8 


0.1 


0.8 


98.3 


459 


0.1 








13.3 


86.6 


335 


0.8 


0.0 


14.6 


84.6 


565 


0.1 





2 


99.5 


0.2 


336 


0.2 


0.0 


63.6 


36.3 


864 


0.0 








99.5 


0.4 


342 


0.0 


0.0 


25.3 


74.7 














351 


0.3 


0.0 


5.6 


94.2 














378 


13.5 


0.6 


77.2 


8.7 














432 


0.0 


0.0 


99.7 


0.3 














528 


0.1 


0.0 


15.1 


84.8 














542 


0.0 


0.0 


3.3 


96.7 














561 


0.1 


0.0 


81.7 


18.3 














622 


0.0 


0.0 


0.7 


99.2 














857 
909 


0.0 
0.1 


0.0 
0.0 


99.6 
0.3 


0.4 
99.6 















TABLE VIL F^u ^ E + Ai_ splitting of the KLT TO soft 
mode for various supercells at fixed a — 3.9879A. Experimen- 
tal results are given for x = 0.05. Frequencies are in cm~^. 



supercell 


X 


uto (E) 


I^TO (Al) 


splitting 


KT40 





(76) 


(76) 





KLT80a 


0.0625 


79 


139 


60 


KLT40 


0.125 


85 


197 


112 


Experiment— 


0.05 


c. 50 


c. 80 


c. 30 
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but the Ai soft-mode component frequency increases by 
only about 0.04 cm"-'^, owing to the small participation 
of Li in this mode. The E-mode frequency increase is an 
order of magnitude smaller. 

Computed, averaged, infrared reflectivity spectra for 
poled KLT40 and KT at the equilibrium lattice constant 
are compared in Fig. 4. The two spectra are generally 
similar, but in KLT40, KT-bands are split by symmetry- 
breaking Li-displacements. 

C. Isolated Li- Ions 

Li-ions can jump between nn [0,0, c?]c wells via 
[1,1, 0]c saddle points. Fig. 5 shows the computed 
Li potential between nn Li-wells; e.g. [0,0, d]c and 
[0, d, 0]c via a minimum energy [0,1, l]c saddle. Li- 
ions were moved along a straight line between the wells. 
The energetics of Li motion, with all non-Li ions fixed in 
the positions that they occupy when Li is in its equilib- 
rium off-center state, are contrasted with the energetics 
for Li motion in conjunction with correlated relaxations 
of non-Li ions. Note that these frozen coordinates dif- 
fer from those in Fig. 2 where all non-Li ions occupy 
positions that are consistent with a center of symmetry. 
Clearly, ionic relaxation promotes Li-hopping by reduc- 
ing potential barriers. Fully relaxed calculations yield 
reasonable agreement with results obtained from Arrhe- 
nius fits to NMR data-. The fixed-ion potential barrier 
is ~103 meV (1190 K). The experimental value, 86 meV 
(1000 K) [3] is lower, but the value 103 meV (1200 K) 
[^^] is close. A similar potential barrier value was also 
obtained in the KLTSOa computation. Note that com- 
puted barrier heights will be systematically larger than 
values obtained by experiment if the latter are reduced 
by Li-tunneling. 



D. Li-Li Pairs 

The energetics of Li-Li nn pair displacement configura- 
tions were studied using the KLTSOa supercell (Fig. 6) 
for all nn-pair configurations listed in Table IVlTll Total 
energy calculations, with all ions relaxed and a = 4 A, 
indicate that the lowest energy path for converting the 
minimum energy z — z configuration into the z — z con- 
figuration is via intermediate microstates; e.g. the en- 
ergy ior z — z ^ X — y is 2627 K (Table I1X|I within 
the experimentally reported range of values, 2100-2800 
j^^nj^ attributed to the rearrangement of Li-Li nn- 
pairs.-'^fi Barrier height computations for two selected low- 
energy paths (see Fig. 7) were performed by converting 
the z — z configuration into z — 'z via intermediate states 
(Fig. 6). Note that the metastable states are separated 
by high barriers the lowest of which, 356 meV (4130 K), is 
significantly larger than the value 2100-2800K^^i^^i^^ ob- 
tained by fitting an Arrhenius expression to experimental 
data; this discrepancy may be explained by Li-tunneling 



TABLE VIII: Fully relaxed Li-Li pair configurations in 
KLT806 (no subscript), KLTSOc (subscript 2), KLTSOd (3) 
and KLTSOe (4). Li-displacement vectors from ideal per- 
ovskite positions are given in A. 



z 


— z 

















X 


— z 


1432 















X 




2570 


1138 













y 


— X 


2627 


1195 


57 











X 


— X 


2990 


1558 


420 


363 









X 


— X 


3170 


1738 


600 


543 


180 







J. 


— z 


3883 


2451 


1313 


1253 


893 


713 





z 


— "z 


6379 


4947 


3809 


3752 


3389 


3209 


2496 



between metastable states that are close in energy. Ionic 
potentials were computed, at fixed lattice parameters, 
in a two-step sequence: 1) A Li-ion in a Li-Li nn pair 
was moved along a straight line between nn equilibrium 
states, while other ionic coordinates (including those of 
the second Li ion) were fully relaxed. 2) The second 
Li-ion was moved along a straight line between nn wells 
while all other ionic coordinates (including those of the 
first Li ion) were fully relaxed. Simultaneous motions of 
both ions in a Li-Li pair were tried, but this path required 
greater energy. Hence, Li dynamics within Li-Li pairs is 
Glauber type^S Significantly, the excitation energy for 
z — z ^ X — X Li-Li nn-pair flipping is close to the en- 
ergy for z — z ^ x — x and z — z — > a; — y flipping. These 



notation 


Xi 


2/1 


Zl 


X2 


J/2 


Z2 


E[eV] 


z - 


2 


0.00 


0.00 


1.43 


0.00 


0.00 


1.19 


0.000 


X - 


- 2 


1.18 


0.00 


0.06 


-0.01 


0.00 


1.31 


0.124 


X - 


- 2 


1.11 


0.00 


-0.02 


0.00 


0.00 


-1.12 


0.221 


V - 


X 


0.00 


1.12 


-0.01 


1.12 


0.00 


0.01 


0.227 


X - 


X 


1.07 


0.00 


0.00 


1.07 


0.00 


0.00 


0.257 


X - 


X 


1.10 


0.00 


-0.01 


-1.10 


0.00 


0.01 


0.279 


2 - 


2 


0.00 


0.00 


-1.09 


0.00 


0.00 


1.09 


0.335 


z - 


2 


0.00 


0.00 


0.76 


0.00 


0.00 


-0.76 


0.550 


X — 


X2 


1.14 


0.01 


-0.01 


1.13 


0.01 


-0.01 


0.141 


X — 


Z2 


1.11 


0.00 


-0.02 


0.03 


0.00 


1.13 


0.213 


X — 


y2 


1.11 


-0.02 


0.01 


0.02 


-1.11 


0.01 


0.216 


z — 


22 


0.00 


0.00 


1.14 


0.00 


-0.00 


-1.16 


0.235 


X — 


J/2 


1.14 


0.03 


-0.01 


0.04 


1.14 


-0.02 


0.258 


z — 


22 


0.00 


-0.03 


1.09 


0.00 


-0.00 


1.09 


0.275 


X — 


X2 


1.10 


-0.02 


-0.00 


-1.18 


-0.05 


0.01 


0.383 


z — 


Z3 


-0.01 


0.00 


1.14 


-0.01 


0.00 


1.14 


0.138 


z — 


X-i 


0.03 


0.01 


1.12 


1.12 


-0.02 


0.03 


0.229 


2 — 


Z3 


0.01 


-0.01 


-1.09 


0.00 


0.00 


1.09 


0.385 


z — 


X4 


-0.02 


0.00 


1.14 


1.11 


0.00 


-0.01 


0.213 


2 — 


24 


0.00 


0.00 


1.09 


0.00 


0.00 


1.08 


0.248 


z — 


24 


0.00 


0.00 


1.10 


0.00 


0.00 


-1.10 


0.257 



TABLE IX: Excitation energies for Li-Li nn-pairs (K). 
notation z — z x — z x — z y — x x ~ x z — z x — x z — z 
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results are inconsistent with the classic dipole-dipole in- 
teraction expression E oc (/2i • /i2) — 3(/Ii • fi2){fi2 ■ ^^12)- 

If the experimental observation of two relaxation 
proccsscE^^-^^ corresponds to flipping of isolated Li-ions 
and Li-Li nn-pairs, then significant concentrations of 
both must be present. Assuming a random distribu- 
tion of Li-ions in KLT, one can calculate the densities 
of isolated Li-ions and Li-Li nn pairs as functions of x] 
x{l — x)^, and 3x^(1 — a;)^", respectively (Fig. 8). 

The average distance between Li-Li nn-pairs is ob- 
viously larger than the average distance between iso- 
lated Li-ions because the concentration of Li-Li nn-pairs 
[xLi-Li) is so much lower (xLi-M ~ 3a;^ for small x). 
Therefore, it is plausible to think of randomly oriented 
but interacting PNC that each contain one Li-Li nn-pair 
as the "seed" plus several isolated Li-ions. Such a sys- 
tem has an interesting phase diagram^^ Because the en- 
ergy required to reorient the polarization axis of a Li-Li 
nn-pair by 90° is large, 2990 K, PNC polar orienta- 
tions resist coarsening that involves reorientation(s) of 
Li-Li nn-pair axes. This constraint on FE-percolationi 
may explain the observation of a dipole-glass like state 
in KLT^ with < a; < 0.02. 

Results for supercells with two Li-ions, separated by 
more than oq, are given in Table Ivml For example, 
KLTSOe has three symmetry-independent configurations 
(Fig. 9): (1) both Li-ions are displaced in the same direc- 
tion, z — Z4; (2) they are displaced in opposite directions; 
z—'zi] (3) they are displaced in perpendicular directions, 
e.g. z — X4. The z — configuration has lowest energy 
at fixed lattice parameter. This may be a consequence of 
accommodating two off-centered Li-ions, each of which 
produces a strong local tetragonal distortion, into a cell 
in which the global strain is constrained to be zero. 

Stachiotti et al^ used a shell model to compute Li- 
Li interactions in KLT out to distance 6 oq. Their re- 
sults also indicate large strain-, and nonlinear polar- 
ization corrections, to the dipole-dipolc interaction en- 
ergy. Experimental diffuse scattering data for KLT*** in- 
dicate strongly correlated planar atomic displacements 
rather than chain-displacements as observed in KNbOs, 
BaTiOa, KTaOg and SrTiOai^-iSi^ 

E. Dynamical Charges and Polarization 

The polarization induced by Li off-centering was stud- 
ied by Berry phase analyses, and dynamical charges 
^aa were determined for all ions in KLT40R,of, fully re- 
laxed KLT40 (Tabic m, and KLT4G with ions in ideal 
perovskite positions f Table IXl)) . Compare these charges 
with Zzz in KT with the same lattice parameter: K: 
1.151; Ta: 8.680; Ox: -1.690; Oy: -1.690; Oz: -6.452. 

The dynamical charge is a symmetric function of dis- 
placements (u) from ideal perovskite positions, so to a 
first approximation, it is quadratic in u: 

z:,{u) ^ z:m ~ (1) 



TABLE X: Dynamical charges Zl^ in the KLT40 supercell. 



Ions 




X 


y 


z 




Zxx 


Zyy 


Ki 


la 


r\ r\r\r\r\ 

0.0000 


0.0000 


0.0014 


1.445 


1.142 


1.142 


K2 


2c 


0.5000 


0.0000 


0.9996 


1.024 


1.180 


1.132 


K3 


2c 


0.5000 


r\ r\f\f\r\ 

0.0000 


0.4963 


1.014 


1.168 


1.122 


K4 


la 


0.0000 


0.0000 


0.4952 


0.999 


1.146 


1.146 


K5 


lb 


0.5000 


0.5000 


0.0074 


1.402 


1.177 


1.177 


Li 


lb 


0.5000 


0.5000 


0.6273 


0.574 


1.106 


1.106 


Tai 


4d 


0.2513 


0.2513 


0.2526 


8.145 


8.745 


8.745 


Ta2 


4d 


0.2494 


0.2494 


0.7499 


8.619 


8.479 


8.479 




AA 


u.zoio 


u.zoio 




— D.UUo 


— l.DoO 


— i.DoO 


02 


4d 


0.2474 


0.2474 


0.4957 


-6.232 


-1.639 


-1.639 


03 


4e 


0.0000 


0.2500 


0.2436 


-1.724 


-6.382 


-1.738 


04 


4e 


0.0000 


0.2533 


0.7457 


-1.661 


-6.358 


-1.619 


05 


4f 


0.5000 


0.2482 


0.2442 


-1.651 


-6.662 


-1.644 


Ob 


4f 


0.5000 


0.2589 


0.7379 


-1.617 


-6.197 


-1.790 



TABLE XI: Dynamical charges Z*^, ionic displacements (in 
A) and total dipole moment (in eA) in the KLT40 supercell. 
Only symmetrically distinct ions are listed: for ionic coordi- 
nates see Table 13 



Ions 


Refer. Ideal perov. Average Displ. 


Tot. dip. 




struct. 


positions 


dyn. ch. 


[A] 


[eA] 


Ki 


1.153 


1.152 


1.250 


0.011 


0.014 


K2 


1.153 


1.151 


1.109 


-0.003 


-0.007 


K3 


1.153 


1.153 


1.107 


-0.029 


-0.064 


K4 


1.153 


1.153 


1.102 


-0.038 


-0.042 


Ks 


1.153 


1.124 


1.217 


0.058 


0.071 


all K 






1.143 


-0.005 


-0.029 


Li 


1.251 


1.251 


1.025 


1.009 


1.035 


Tai 


8.713 


8.710 


8.522 


0.021 


0.706 


Ta2 


8.714 


8.712 


8.681 


-0.001 


-0.035 


all Ta 






8.602 


0.010 


0.670 


Oi 


-6.450 


-6.472 


-6.315 


-0.051 


1.293 


O2 


-6.529 


-6.504 


-6.413 


-0.034 


0.867 


all 0|| 






-6.364 


-0.042 


2.160 


O3 


-1.703 


-1.700 


-1.708 


-0.051 


0.346 


O4 


-1.700 


-1.700 


-1.687 


-0.034 


0.231 


O5 


-1.687 


-1.685 


-1.673 


-0.046 


0.306 


Oe 


-1.687 


-1.685 


-1.662 


-0.096 


0.637 


all Ox 






-1.683 


-0.057 


1.521 


all 










3.681 


Total 


5.357 



Exact value 5.153 



The average dynamical charge, which is required for 
calculating the polarization, is 

ui 

{z;j = / z:,{u)du = 

ui-uo J 

Uq 
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= Z:M + Ic{uI + uI + u,uo) (2) 
In particular, if uo = then 



{ZL)^^izL{^i) + 2z:M) (3) 

Average values of the dynamical charges, and individ- 
ual ion contributions to the total KLT40 supercell dipole 
moment are listed in Table Ixll For comparison, the ex- 
plicit value of the total dipole moment, obtained by Berry 
phase analysis, is also shown. 

The largest contribution to the total dipole moment is 
from O ions. The Li-displacement contribution is roughly 
equal to that of Ta. The total dipole moment is ^5 times 
larger than the dipole moment due to the Li ion itself, 
which indicates that the Li-dipole moment is strongly 
enhanced by structural relaxation. 

A similar enhancement was also obtained in a shell 
model calculationj^ but in that study the main cause 
of enhancement was a large Ta-displacement. Here, en- 
hancement of the Li-dipole moment is primarily caused 
by coupling between Li- and 0-displacements. 



TABLE XII: The KLTSOb supercell dipole moment (in eA) 
for different Li-Li nn-pair configurations 



notation 


direct. 


K 


Li 


Ta 


Ox 


Oil 


Total 


z 


— z 


z 


0.138 


2.688 


4.636 


2.215 


3.850 


13.527 


X 


— z 


z 


-0.170 


-1.168 


-3.702 


-1.432 


-2.524 


-8.996 






X 


0.335 


1.138 


4.915 


1.085 


1.888 


9.361 


X 


— X 


X 


0.393 


2.184 


5.681 


2.079 


3.417 


13.754 


X 


— z 


z 


0.152 


1.408 


3.830 


1.187 


2.202 


8.779 






X 


0.512 


1.193 


6.315 


0.580 


1.008 


9.608 


V 


— X 


y 


0.393 


1.147 


5.451 


0.867 


1.480 


9.338 






X 


0.393 


1.147 


5.451 


0.867 


1.480 


9.338 



Total dipole moments for a KLT806 supercell with dif- 
ferent Li-Li nn-pair configurations are listed in Table lXIll 
They were calculated with the average dynamical charges 
listed in Table DCII and ionic displacements relative to ideal 
perovskite. Individual ionic contributions are also shown. 
Oxygen ion contributions are separated into 0|| and O^, 
with the polarization direction parallel and perpendicu- 
lar to the Ta-O-Ta bond, respectively. These results indi- 
cate that the total dipole moment induced by Li-Li pairs 
is larger than that from two isolated Li-ions. 



F. Coupling Between Strain and Polarization 

The influence of applied stress on dipole moment and 
energy was studied by applying various tetragonal dis- 
tortions to the KLT40 supercell with a = 3.95594 A, and 
0.98a < c < 1.02a. Total energy calculations show that, 
for c > a, configurations with [0, 0, l]c-displaced Li-ions 



are lower in energy than those with [1, 0, 0]c-displaced Li. 
The difference in energies is 

- 5E,^^ [meV] = 1.649e3 -f 105. 7e^ -f 3430.e^ (4) 

where 63 = (c/a — 1). By using this energy difference, 
the occupation probability for Li-displacements in the c- 
direction can be estimated as 

« 1/6 - 5E,J%kBT (5) 

Hence, tetragonal strain increases the population of Li- 
displacements parallel to the c-axis. The energy decrease 
(per primitive cell) connected with redirected Li-ions is, 

< SElizzzz >= x5Ezz{W-z-z + Wzz ~ 1/3) 

« -2x5ElJQkBT (6) 

which at small 63, is proportional to e§. Thus, the en- 
ergy decreases by local Li-redirections along c and this 
decreases the corresponding elastic constant proportion- 
ally to X. 

Li-redirection also increases the square of supercell 
dipole moment dz along the c-axis 

5qzzzz = d^dz{wzz + w-z-z - 1/3) 

^ 2dlSEzz/9kBT. (7) 

By comparison with Sqzzzz is proportional to strain 
for small 63. Because Li- redirection is coupled to strain, 
one expects a relaxation contribution to the electrostric- 
tion constants and acoustic response, as observed by Pat- 
tnaik and Toulouse fiii^]. 

In addition to promoting Li-redirection, tetragonal dis- 
tortions change the TO soft mode frequency. Using 
the average dynamical charges and computed ionic dis- 
placements, one can calculate supercell dipole moments, 
dz and d^ in eA, which correspond to Li-ions displaced 
in the [0,0, l]c and [1,0, 0]c directions respectively, for 
a < c < 1.2a: 

4 = 5.15 -M67.e3 + 2987.e^ 

= 5.15-40.63 -1- 585.6^ (8) 

These results show that the supercell dipole moment 
is greatly enhanced by an increase in c/a. Partial con- 
tributions to the dipole moment that are associated with 
Ta- and 0-displacements are most sensitive to changes 
in c; which is consistent with an electrostriction induced 
reduction in the soft-mode frequency. 
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G. Toward an Effective Hamiltonian 

The picture of KLT in which Li-dipoles, a soft-mode, 
and PNC all interact, lends itself to a first-principles 
based effective Hamiltonian, Heff, treatment. Full for- 
mulation of Hf,f f is beyond the scope of this work, but 
some essential intermediate results have been derived. 
For simplicity, consider KLT40, with the same ionic dis- 
placements in each unit cell. 

To derive He// for a FE, one selects a local basis 
for the instabilities that cause the FE transition, and 
includes interactions between sites, between local distor- 
tion and strain, etci^. In KLT, it is clear that instabil- 
ities centered on Li are responsible for the phase tran- 
sition, and that Li-centered distortion variables must be 
included in H^f f . However, the coupling of Li-dipole mo- 
ments to the KT soft-mode is sufficiently important that 
one must also include a local variable for the KT soft- 
mode. By analogy with -ffe// for KNbOg^, this local 
variable is centered on the Ta sites. 

Several eigenvectors are relevant: 

• vi^a ■ the dynamical matrix eigenvector for the KT 
soft-mode, polarized along a. 

• W2,a '■ the dynamical matrix eigenvector for the Li- 
instability of KLT40, polarized along a. (Dynami- 
cal matrix eigenvectors are normalized for 40-atom 
cells) . 

• VTa,a = vi^a and ULi,Q = V2,a, orthogonalizcd to 
Vi^cn and renormalized. 

• dLi,a and dTa,a, the dimensionless displacement 
eigenvectors obtained by dividing the elements of 
VTa,a and vi^i aj respectively, by {m/nio), where 
(arbitrarily) mass toq = 1 amu. 

Consider displacement patterns in KLT40 in which 
displacement amplitudes corresponding to VTa,a and 
VLi.a are Tq, and Aq, respectively. The energy is mini- 
mized with respect to all other modes. The strain tensor 
is {ei}, i = 1,6, in Voigt notation. Calculated results 
closely fit the expansion: 

U = Uo 

- 0.055566|Ap + 0.004160|A|* 

+ 0.016303(A2a^ + A^A^ + A^A^) 

- 0.000077|A|^ + 0.000001|A|^ 
4- 0.024760|tP + 0.008677|t|* 

- 0.011288(rjr2 -f rjr^ -|- r^rf ) -|- 0.000684|t|'5 

- 0.025236A • f + 0.001210(A^t^ + \\Ty + A^r^) 

- 0.000370(A^t^(A^ + A^) cyclic perm.) 
817.342(6? + e2-Fe2)_^ 

+ 210.580(61624-6163-1-6263) 



-f 732.803(6^ 6^ + 6^) 

~ 0.262724(6iA2 +62A2 +63A2) 

- 0.209816(6i(A2 + a2) + 62(A2 + A^) 
+ e3(A2+A2)) 

- 5.30307(64Aj^Az -I- e'^\x'>^z + e^X^Xy) 

- 3.24098(6itJ + 62t2 -f e^Tl) 

+ 0.715628(6i(Ty2 + tI) + e^ijl + t^) + 
+ e3(rJ+T2)) 

+ 1.38297(64ryr^ -I- c^t^t^ + eeT^Ty), 

(9) 

where U is in eV and A and r in A. The con- 
nection between TO polarization and r is, PtOz = 
3.766820x^/1^ eA where V is the unit cell volume. 
The connection between the Li-dipolc moment and A is, 
= 0.876855A^ 6A. 

The potential for A alone (~ Li off-centering) has 6 
wells in [0,0, l]c directions. The potential for f alone 
(KT soft-mode) has a small positive harmonic coefficient 
and is therefore highly polarizable. Negative bilinear cou- 
pling between A and r indicates that the soft-mode en- 
hances the total polarization from Li off-centering. Posi- 
tive coefficients of |r|^ and |t|^ imply that the soft-mode 
frequency hardens as the structure distorts. Higher-order 
coupling between A and r is included because it signif- 
icantly improves the fit. 

One can reduce the number of variables by excluding 
strain. This results in a decrease of the nonlincarity con- 
stant which multiplies the fourth power of polarization, 
and changes the sign of the angle dependent term for the 
soft-mode. 

Eq. Inland the mode effective charges can be used to cal- 
culate the soft-mode and Li-dipole contributions to the 
susceptibility. The soft-mode contribution is 110 for the 
reference state, much less than might be expected con- 
sidering that KT is an incipient ferroelectric. The LDA 
overestimates the soft mode frequency of KT, which leads 
to too large a value for the jrp coefficient of Eq. Inland 
too low a permittivity. For fully relaxed KLT40, the 
model predicts that the contribution of the soft-mode to 
the susceptibility along the polar axis is about 55 due to 
mode hardening. The Li-dipole contribution is insignifi- 
cant (order 1). This, however, is only the contribution to 
the static susceptibility from Li vibrating within a sin- 
gle well. As discussed in Section HVI Li-hopping gives a 
much greater, frequency-dependent, contribution to the 
susceptibility. 

From Eq. |51 and the mode effective charges, the 
coupling between TO phonons and Li-dipoles is 
0.064:fizPTOz/iVeo where £0 is the permittivity of free 
space. This coupling enhances dielectric permittivity and 
shifts the Curie temperature to higher values (c.f. Sec 
II V|) . This Lorentz factor, 0.064, is smaller than that 
found''* in a point polarizable ion model, (0.1 - 0.2). The 
Lorenz factor decreases quadratically with Li displace- 
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ment and equals 0.037 at the equilibrium Li position. 

To complete H^ff requires: (1) Formulating lo- 
cal representations of A and r e.g. lattice Wannier 
functions^); (2) Quantifying distance dependencies of 
the interactionsSiiSi. 



H. Band Structure Calculations 

Strong n-type photocurrent is observed in KLT at 
T < 80K^^, and it has been attributed to shallow 
Oe2p hole states at the bottom of the valence band which 
are caused by large Oe-displacements (notation of Table 
Ivji . Holes are trapped by these levels which prevent re- 
combination with electrons that were promoted to the 
conduction band by light absorption. 

To study shallow forbidden gap states in KLT, band- 
structure calculations were performed for KLT40 and the 
corresponding ideal KTaOs supercell, in the T ~ X di- 
rection. Results for the top of valence and bottom of 
conduction bands are shown in Fig. 10. The top of the 
valence band, and the bottom of the conduction band, 
are split and/or distorted in KLT40 relative to KT. 

VASP results confirm large Oe-displacementS)^ and 
demonstrate some splitting at the top of the valence 
band, and some distortion at the bottom of the conduc- 
tion band. This occurs because Li off-centering breaks 
Pm3m symmetry creating symmetrically different 0-ions 
with different self-consistent electrostatic potentials (cf. 
Rcf.'"'^). The VASP results, however, do not confirm the 
presence of OgSp states at the top of the valence band. 



TABLE XIII: Electrostatic potentials (in eV) on the ions of 
the KLT40 supercell in the equilibrium structure (Equil) , ref- 
erence structure (Refer), and in pure KT (KT). Only sym- 
metrically distinct ions are listed (atomic positions are listed 
in Table El 



Ion 


Equil. 


Refer. 


KT 


Ion 


Equil. 


Refer. 


KT 


Ki 


10.15 


10.18 


10.33 


K2 


10.13 


10.15 


10.33 


K3 


10.18 


10.23 


10.33 


K4 


10.15 


10.15 


10.33 


Ks 


10.14 


10.23 


10.33 


Li 


43.90 


42.68 




Tai 


-1.21 


-1.22 


-1.30 


Taa 


-1.29 


-1.22 


-1.30 


Oi 


-56.27 


-56.26 


-56.07 


O2 


-56.30 


-56.30 


-56.07 


O3 


-56.19 


-56.26 


-56.07 


O4 


-56.46 


-56.26 


-56.07 


O5 


-56.28 


-56.31 


-56.07 


Os 


-56.43 


-56.31 


-56.07 



Computing the electrostatic potentials on different ions 
in KLT40 (Table IXlTHl reveals that the 062p states are 
significantly deeper in fully relaxed KLT40 than in KT 
or KLT40Rof- The same result was obtained from pro- 
jections of band states onto Oe2p states: 062j>-states in 
KLT40 lie below the top of the valence band. Similar 
results were obtained for Li-Li nn pairs in KLT80, where 
the tetragonal distortion from Li also causes valence and 
conduction band splitting. 



IV. PHENOMENOLOGY 
A. Li-pair relaxation time 

Consider KLT as a distribution of PNC (containing 
one or more Li-ions) in a KT lattice. FP calculations 
on isolated Li-ions and on Li-Li pairs have shown that 
each PNC may have a set of states defined by the Li- 
displacement directions. For a given PNC, group theory 
can be used to enumerate distinct relaxation processes 
and identify which of them couple to an electric field. 
The kinetic equations for state occupation probabilities 
can be solved to estimate the relaxation times. If PNC 
do not interact, then there are distinct relaxation times 
for each cluster type, and the relaxational contribution 
to the dielectric function is the sum of contributions from 
each cluster. Interactions between clusters will broaden 
an otherwise sharp distribution of relaxation times. 

For isolated Li-ions, there are six states, and symmetry 



analysis of relaxation processes yields A 
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En 



Flu- 



The Aig "process" is the trivial steady state equilibrium 
solution. Only the Fiu process couples with an electric 
field and affects dielectric dispersion. The Eg process 
couples with stress, but not electric field. Let k be the 
hopping rate between wells separated by 90° and ne- 
glect hopping between wells separated by 180°. Then 
relaxation time = Ak for the Fiu and r"^ = 6fc for 
the Eg processS,. 

For Li-Li nn-pairs, there are 36 states, many of which 
are metastable (Table HXjl . Symmetry analysis identifies 
seven relaxation processes that contribute to the dielec- 
tric function. Six of these involve substantial probability 
redistribution among metastable states whose occupan- 
cies are negligible. Thus, only one relaxation process 
contributes significantly to the dielectric function; the 
one whose net effect is z — z ^ z — z. 

The full set of kinetic equations give a 36 x 36 matrix. 
For simplicity, consider only the states that are associ- 
ated with the minimum energy path in Fig. 7: 



dwg 
dt 
dw\j 

'~dr 

dwc 

'~dr 

dwd 

'~dr 

dwe 



■ kabWa - haWb 
{ha + hc)Wb - kabWa - k^bWc 
2kcbWc - kbcWf, - hcWd 
{kbc + kba)Wd - kchWc - kabWe 
kabWe " kbaWd (10) 



where Wa is occupation probability for well a. The re- 
laxation time of the dominant process is: 



{kab + ha + he) 12 - \/ a?- - kabk 



ba ; 



(11) 
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Here a = {kat + ha + kbc)/'2- 

In the realistic limit kab ~ fcbc ^ fccb ~ ^6a: r^^ ~ 
kbckab/kba- This rate corresponds to overcoming the 
principle Li-Li nn-pair barrier {z — z — > a; — in the 
z — z ^z — z process. 

B. Dielectric-permittivity frequency-dependence 

If one excludes strains from a mean field Hamiltonian 
then there are two variables: one for a soft-mode; the 
other for Li-dipoles. Both contribute to the dielectric re- 
sponse, but the former is essentially dispersionless far be- 
low the soft-mode frequency. The Li-dipole term however 
has relaxational character. Neglecting the damping of Li- 
oscillations in Li wells, but accounting for jumps between 
nn Li-wells, the Li-dipole associated with z-polarization 
of an isolated Li-ion in a cubic environment is: 

= ^ {yj.. - w^,) (12) 

where V is the primitive unit cell volume. The Fiu re- 
laxation process contributes to polarization fiuctuations, 
and the polarization takes the form^ 

= eoF{T)E, (13) 

where 



Li-dipoles and the TO soft mode; Tc increases and fre- 
quency dispersion is enhanced. These results were ob- 
tained by considering all hopping paths between the six 
[0, 0, djc-type wells, not by invoking an Ising model. 

The dominant relaxation process for Li-Li nn pairs 
leads to an additional relaxational contribution to F{T), 
proportional to Li-Li nn-pair concentration, and with 
T^^ = kbckab/kba (cf. results in Refs»i2iS). The concen- 
tration of Li-Li nn pairs is low, but their dipole moments 
are large (Table IXII() , and xli is proportional to the 
dipole moment squared l|14|) . 

The temperature dependence of the kinetic coefficient 
k is common, and it is mainly controlled by the barrier 
height U. The temperature dependence of the relaxation 
time is 

oo 

T =Ax{T) J e(''+^y")^''-^dy = 

— OO 

^Ax(T)y^e^(^)/'=-^ = 

= ro(T)e^(^)/^-^, (17) 

where x is the dielectric permittivity of the media in 
which the PNC are embedded and a is a constant deter- 
mined by the force constants at the saddle point. 

V. CONCLUSIONS 



F{T) 



1 



1 



Veo eksTAk- 



lUJ 



(14) 



Here: Ez = + XPtoz is the local field, A = Aq — 
3Ai/i^ a coupling constant, Ao and Ai are the effective 
Hamiltonian coefficients of ^fJ-zProz and ^i^Proz re- 
spectively; Ez is an external field; Ptoz is the polar- 
ization connected with the soft-mode; and r = 4fc is 
the relaxation time for the Fiu relaxation process. The 
dielectric susceptibility connected with Li jumps can be 
now written aa^i^SiS 



XLi 



F{T) 



1 - A2F(r)eoM(T) 



(15) 



where 



A{T) = a + 3/3P^oz + ^iPtOz 



(16) 



Here the nonlinearity coefficient (3 is renormalized af- 
ter the elimination of strains; a is renormalized due to 
fields produced by Li. The positive correlation between 
a and increased x implies hardening of the soft mode 
with increasing Li concentration. Dielectric susceptibil- 
ity is enhanced at low frequencies by coupling between 



Calculations on a 40-atom KLT40 supercell that con- 
tains one Li-ion indicate that there is only one instability 
connected with Li off-centering. This instability is associ- 
ated with large Li-displacement in the [0, 0, l]c direction. 
Li potential wells are separated by potential barriers with 
[1, 1, 0]c saddles. Calculated [1, 1, 0]c potential barriers 
and [0, 0, d]c Li-displacements, in fully relaxed KLT, are 
in good agreement with values obtained by fitting a mean 
field expansion to experimental datciii*i^. 

All ions in the KLT40 supercell are displaced from ideal 
perovskite positions, but the 0-ions closest to Li exhibit 
the largest sympathetic displacements, towards Li. Cal- 
culated vibrational frequencies for fully relaxed KLT40 
indicate that Ai- and i?-symmetry modes emanate from 
an Flu TO mode similar to the soft mode of pure KT. 
The average frequency of these modes is higher than that 
of the TO soft mode, which reflects hardening of the soft- 
mode induced by substituting Li for K. The intrawell Li 
vibrational frequency is about 300 cm^^. The soft-mode 
and Li-polarization are coupled. 

Berry phase analyses yielded total dipole moments as- 
sociated with Li off-centering that are ^5 times larger 
than the Li-dipole moment itself. This enhancement is 
primarily caused by 0-relaxation around the off-centered 
Li-ion. 

Fully relaxed KLT40 is tetragonal with c/a > 1, and 
polarization along the c-axis. Polarization increases with 
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increasing c/a; uniaxial strain. grant # 01-02-16029 and thanks NIST for hospitality. 

First-principle calculations of Li-Li nn-pair excitations 
in a KLT80 supercell strongly support the conclusion of 
Dousenau et alm^ that the "180°" relaxation process is 
caused hy z — z^'z — 'z Li-Li nn-pair flipping. 
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Captions 

Fig. 1. KLT configurations compared on a common 
320-atom supercell. The grid connects perovskite A-sites. 
Empty sites represent K; black circles represent Li. When 
there is one Li per primitive cell, it is arbitrarily displaced 
in the z direction; when there are two Li per primitive 
cell, one is displaced in the z direction and the other in 
the X direction. 



Fig. 2. Structural distortions around Li in KLT40. 

Fig. 3. Energy as a function of Li-displacement from 
an ideal A-site position in KLT40: with all ions' coordi- 
nates frozen in the centrosymmetric positions except for 
Li (a); with full relaxation of all ions' coordinates except 
for Li (b). 

Fig. 4. Comparison of the reflectivity spectra com- 
puted for KT and KLT40 (the component along the polar 
axis {Rz)^ perpendicular to it (Rx), and averaged spec- 
trum (Rav)) at identical lattice parameters ao — 3.956 A. 
The electronic permittivity £oo was set to 5.15 in each 
case, the experimental value for KT^ The damping con- 
stant was set to 20 cm~^ for all modes in KT and 40 
cm^i for KLT40. 

Fig. 5. Li potential barriers in KLT40 in the direction 
between nearest neighbor Li-wells computed with all ions 
except Li fully relaxed (relaxed) , and with all ions except 
Li frozen at their coordinates corresponding to Li at its 
equilibrium position (frozen). 

Fig. 6. A-site configurations in KLT805. Undecorated 
sites are occupied by K. 

Fig. 7. The energetics of Li-Li nearest neighbor pair 
reorientation, z — z — > ^ — z, including the intermediate 
states depicted in Fig. 6. 

Fig. 8. Densities per A-site of isolated Li-ions and Li- 
Li nearest-neighbor pairs in Ki-^Li^^TaOa, as functions 
of X. 

Fig. 9. Fourth nearest neighbor Li-Li configurations: 
a) {z — Z4); b) (z — Z4), and c) (z — 2:4). 

Fig. 10. Comparison of the electronic band structures 
in KT40 (a), KLT40 (b), and in KLT80 with a z - z Li 
nn-pair (c). The valence states are below zero energy. 
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FIG. 1: Prosandeev, Cockayne, Burton, PRB 




FIG. 2: Prosandeev, Cockayne, Burton, PRB 




FIG. 4: Prosandeev, Cockayne, Burton, PRB 
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FIG. 7: Prosandeev, Cockayne, Burton, PRB 




FIG. 8: Prosandeev, Cockayne, Burton, PRB 



